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ABSTRACT 

Wc study numerically the relativistic Bondi-Hoyle accretion of an ideal gas onto a 
Kerr fixed background space-time on the equatorial plane with s-lab symmetry. We 
use both Kerr-Schild (KS) and Boyer-Lindquist (BL) coordinates. We particularly 
focus on the study of the flip-flop motion of the shock cone formed when the gas is 
injected at supersonic speed. The development of the flip-flop instability of the shock 
cone in the relativistic regime was reported recently for the first time. We reproduce 
the flip-flop behaviour found in the past when BL coordinates are used, and perform 
similar numerical experiments using horizon penetrating KS coordinates. We flnd that 
when using KS coordinates the shock cone oscillates, however such oscillations are not 
of the flip-flop type and their amplitude decrease with resolution. 

Key words: Accretion, Accretion discs - shock waves - hydrodynamics - black hole 
physics 



1 INTRODUCTION 

Bondi-Hoyle accretion is the evolution of a homogeneously 
distributed gas moving uniformly toward a central compact 
object ( |Bondi fc Hoyle et al. 1944] ). Depending of the rela- 
tion existing between the velocity of the gas and the sound 
speed a shock cone can be formed. Specifically, when the 
fluid is supersonic, a shock cone appears, otherwise it does 
not. 

When the compact object is rotating, an interest- 
ing phenomenon, called flip-flop instability can appear 
(|Matsuda et al. 1987|l . This effect consists in the oscillation 
of the shock cone orientation. Various studies of this effect 
have been carried out under different conditions and for dif- 
ferent purposes. For example, in Newtonian gravity, consid- 
ering the density and velocity field to be non-uniform, they 
found that the shock cone might achieve a disk-like topology 
( |Fryxell fc Taam et al. T988| |Taam fc Fryxell et al. 19881 ); 
also the dependence of the appearance of the flip-flop be- 
haviour in terms of the size of the object has been studied 
(|Sawada et al. 1989: Mats uda et al. 19911 ILivio et al. 19911 
IBenensohn et al. 19"97| . 

The astrophysical relevance of the flip-flop of the shock 
cone is associated for example with the fluctuations observed 
in the X-ray emission of pulsing systems such as Vela-xl and 
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EXO 2030+375, that are not associated to the orbital period 
of the binary system ( [Taam fc Fryxell et al. 1988[ ). 

The flip-flop effect can also happen within a relativistic 
scenario. Particularly interesting is the case when the com- 
pact object is a black hole. This case was first reported in 
()Donmez et al. 20111) and we consider this is a very impor- 
tant subject to be studied in detail, because the oscillations 
(of different nature) of a gas around a black hole may po- 
tentially explain the oscillatory behaviour of high energy 
sources, like quasi periodical oscillating sources (QPOs) or 
bursts when the shock cone is unstable. Moreover, we con- 
sider the treatment involving the curved background space- 
time to be a very important step, which uses a truly black 
hole space-time background instead of the commonly used 
pseudo-Newtonian Paczyhski-Wiita type potentials. 

We are in the position of analyzing the case where the 
gas is truly allowed to enter the black hole's event hori- 
zon, this is what we do in this paper. In the relativistic 
Bondi-Hoyle accretion, the flip-flop instability in the equa- 
torial plane was briefly reported by (|Donmez et al. 20111) . 
However, the coordinates used to describe the rotating Kerr 
black hole are not appropriate, because the Boyer-Lindquist 
coordinates are singular at the event horizon, and as a con- 
sequence it is necessary -as in the Newtonian cases- to apply 
boundary conditions at an inner time-like artificial bound- 
ary outside and near the horizon. 

In this work we study the Bondi-Hoyle accretion process 
using the Kerr-Schild penetrating coordinates in the equato- 
rial plane, in these coordinates the fluid naturally falls into 
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the black hole. We also study the case of using the singular 
Boyer-Lindquist coordinates in order to reproduce the flip- 
flop effect and measure its amplitude and phase. We repro- 
duce the morphological behaviour in (|Donmez et al. 201fp . 
and measure the amplitude of the flip-flop oscillations using 
Boyer-Lindquist coordinates. However, when using KS co- 
ordinates, there are some oscillations with a frequency that 
is not comparable to that observed when using BL coordi- 
nates associated to the flip-flop effect; we observe that these 
oscillations are due to numerical artifacts because their am- 
plitude converges to zero when resolution is increased . 

In fact, there is a precedent within the Newtonian 
regime where it has been found that shock cone instabili- 
ties can be associated with the implementation of boundary 
conditions at the surface of the accretor, and may also in- 
fluence the attached or detached nature of the shock to the 
compact object or numerical parameters in the calculations 
( |Foglizzo et al. 2005p . 

This paper is organized as follows. In Sec. [2] we describe 
the system of the relativistic hydrodynamics equations and 
the metric functions for each coordinate system, in Sec. [3] we 
describe the numerical methods used for evolving the equa- 
tions, initial wind conflgurations and boundary conditions. 
In Sec. Uwe have a comparison of the flip-flop behaviour be- 
tween Boyer-Lindquist and Kerr-Schild coordinates system 
and in Sec. [5] we present a discussion of our results and some 
conclusions. 



ds^ = -(a - I3il3'')dr + 2Pidx'^dt + ■ji^dx^dx', (3) 

where /?' is the shift vector, a the lapse function and "/ij axe 
the components of the spatial 3-metric. 

The background space-time corresponds to a rotating 
black hole, however, since the main aim of the paper is the 
study of the influence of a time-like internal artificial bound- 
ary used in previous studies we consider the case of the black 
hole being described with BL coordinates which is the usual 
approach, and the case of using KS coordinates, which have 
the property of penetrating the black hole event horizon. In 
both cases we describe the space-time using spherical coor- 
dinates for the space x'^ = (t, r, 0, (^). 

In the first case, for BL coordinates the metric functions 
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2 RELATIVISTIC HYDRODYNAMIC 
EQUATIONS 

We model the wind using the relativistic Euler equations for 
a perfect fluid described by the stress-energy tensor 



(1) 



where po is the rest mass density of a fluid element, u'^ the 
4- velocity of the fluid elements, h = 1 -\- e -\-p/ po is the spe- 
cific internal enthalpy, with e the rest frame specific internal 
energy and p the pressure. The space-time background is de- 
scribed by the metric , and we use units where G = c = f 
from now on. Additionally, we assume the fluid obeys an 
ideal gas equation of state given by 



p = poe(r-f), 



(2) 



where V is the adiabatic index, or the ratio between specific 
heats. Relativistic Euler's equations are derived from the 
local conservation of the rest mass density and the local 
conservation of the stress-energy tensor 
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and M, a are the mass and the angular momentum of the 
black hole respectively. 

On the other hand, in KS coordinates the metric func- 
tions are 
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Vf.ipou") = 0, 

v„(r'"') = 0, 

where is the covariant derivative consistent with 

Qfj,^ (|Misner et al. 1973)) . It is convenient to cast 
these equations in a fiux balance law fashion, as de- 
scribed in (|Marti et al. 19911 [Banyouls F. et al. IQQT) 
IFont et al. 2000|) . The only ingredient left is the space-time 
background written within the 3+1 decomposition of the 
space-time formalism. Such metric is written in general as 



At this point we assume the system is cylindrically sym- 
metric, that is, we only study the evolution of the gas on the 
equatorial plane and set 6 — 7r/2, which does not mean -as 
usually assumed- this is a sort of an infinitesimally thin mor- 
phology approximation. 

With the space-time background metric known, rel- 
ativistic Euler equations on a curved fixed background 
space-time in our coordinates and for the two cases 
considered read (|Marti et al. 19911 [Banyouls F. et al. 1997[ 
IFont et al. lOOOll 
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dtU + drF'' + d^.F'^ = S - ^^F'' 



V7 



(7) 



where 7 — det{'yij) is the determinant of the spatial metric, 
U is the vector whose entries are conservative variables that 
depend on the original primitive variables po,v^,p,e. The 
vector contains the fluxes along the spatial directions and 
S is the source vector. All these ingredients are specifically 
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where are the components of the velocity of the gas 
measured by an Eulerian observer, which are related to 
the spatial components of the 4-velocity of the fluid el- 
ements by v'' = /W + P' /a, W is the Lorentz factor 
W = — "fijV^v^; are the Christoffel symbols ob- 

tained from the space-time metric. 

Finally, the system of equations (|7I8I9I10I1H) is a set of 
four equations for the variables po, v'~ , e. Thus, in order 
to close this system of equations, we consider the fluid obeys 
the equation of state ([2]). 



3 NUMERICAL METHODS 



We evolve numerically the system of equations (|7l8l9ll0llip 
as a 2D initial value problem using a finite volume ap- 
proximation on the equatorial plane. Specifically we use the 
well known high resolution shock capturing methods, with 
the HLLE approximate Riemann solver (|Harten et al. 19831 
lEinfeldt. 1988|) . and the minmod linear piecewise recon- 
structor. For the update in time we use the method of lines 
with a second order Runge-Kutta TVD integrator along the 
time direction. 



3.1 Initial data 

As initial data, we consider a homogeneous wind, that uni- 
formly fills the whole domain, moving on the equatorial 



plane along the x direction with constant density and pres- 
sure. The initial velocity field can be expressed in terms 
of the asymptotic velocity follows 
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and the relation v'^ = ViV^ — is satisfied. Specifically, in 
the case of the BL coordinates the above expressions for the 
velocity field reduce to 

(18) 
(19) 

In order to choose the initial pressure profile, we introduce 
the asymptotical speed of sound Csoo ■ Once we fix the value 
of Csoo and assume the density to be a constant po — Pini , 
the pressure can be found from the following expression 




Pini 



3 Pini (r — 1) 



r(r- 



(20) 



As we can see, in order to avoid negative and zero pressures, 
the condition Csoo < Vr — 1 has to be satisfied. Finally, the 
internal specific energy is reconstructed from the equation 
of state (0 . 

We find useful to introduce the asymptotic Mach num- 
ber in terms of Vao and Csoo as Moo = This quantity 
determines whether the flow is supersonic or subsonic. When 
Mao is bigger than 1 the flow is supersonic, otherwise the 
flow is subsonic. In this paper, we are only interested in the 
supersonic case as we want to study the properties of the 
shock cone formed, in particular the possibility of flip-flop 
oscillations. 

Other parameters we have to take into account are the 
black hole mass Af and its angular momentum a. In the case 
of the black hole mass we chose units in which M = 1. On 
the other hand, we chose various representative cases for the 
rotation parameter a. In Table 1, we summarize the set of 
initial parameters we consider in our analysis. 



3.2 Evolution 

The numerical domain is G [0, 27r) and r G [rexc,rmax]- 
The election of both the interior boundary at r = r^xc and 
the exterior one at r = Vmax deserves a careful analysis. 
First of all, the interior boundary may be space-like or time- 
like, depending on whether it lies inside the black hole event 
horizon or outside respectively. In most of previous analyses 
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Table 1. In this table we present the models studied in this 
paper where the models are labeled with BL and KS. In the ini- 
tial fluid configuration we have six free parameters, in particular 
we choose three values of the Mach number for supersonic fluid 
configurations. We choose two values of the angular momentum 
of the black hole and two different values of the adiabatic index. 
Notice that the numerical excision boundary is always inside the 
event horizon when using KS coordinates and outside of it when 
using BL coordinates, and the outer boundary of the domain in 
the r coordinate depends on the accretion radius race- The radii 
race, rcxc and r^ax are in units of M. In all the models wo set 
the initial rest mass density to po = 10^® and the asymptotic 
value of the sound speed c^ao = 0.1. 



cal technique used in numerical relativity called black hole 
excision (|Seidel fc Suen 1992|) . It is clear that even though 
the spatial slices penetrate the horizon we cannot include 
the singularity as a part of the numerical domain, and we 
choose < Tcxc < Teh- This implies indirectly that the 
election of r^xc will depend on the rotation parameter of the 
black hole a that defines the event horizon radius teh- 

The possibility to choose r^xc inside the event horizon 
is the core of this paper, because we want to investigate 
whether or not the flip-flop instability of the shock cone in 
the relativistic Bondi-Hoyle accretion is due to the use of 
inappropriate coordinates. 

With respect to the external boundary, the election 
of Tmax IS also restricted. The formation of a shock cone 
strongly depends on the condition that on the one hand the 
gas has to move at supersonic speed and on the other, not 
all the material has to be accreted by the black hole, which 
can happen by simply considering race < rmax, where race 
is the radius of accretion within which all the material is 
captured by the black hole and is defined in relativity by 
(|Petrich et al. 1989|) 

= 75 I ,,2 • (21) 

The boundary conditions used at this boundary are out- 
flow in the downstream region and we inject gas with 
the initial density and velocity in the upstream region 
()Matsuda et al. 19871 IFont fc Ibanez et al. 1998a|> . 

The numerical grid is uniformly spaced in the r and (j) 
directions. In particular we use the base resolution of Ar = 
0.033 in r direction and Atj) = 0.014 in the direction for 
the evolution of all the models. We use a constant time step 
given by At = Cmin{dr, d(fi), where C is a flxed and constant 
Courant factor estimated empirically C = 0.25. Finally, in 
order to avoid the divergences in our variables along the 
evolution, due to the deflnition of the speciflc enthalpy, we 
introduce an artiflcial atmosphere, that is, we impose the 
rest mass density to be no smaller that the minimum value 



such boundary is time-like because the coordinates used are 
Boyer-Lindquist, which are singular at the event horizon, 
and as a way to avoid dealing with the coordinate singu- 
larity there, rexe is chosen simply to be outside the event 
horizon; the price to pay with such an election is that this 
internal boundary requires the implementation of boundary 
conditions for the fluid variables in a region where the grav- 
itational fleld is strong and the gradients of the variables 
involved are big, which eventually implies the propagation 
of numerical errors. 

Our approach involves the use of KS horizon penetrat- 
ing coordinates for a rotating black hole. What can be done 
is to define rexe < rsH, with rEH the radius of the event 
horizon at the equatorial plane, and defines a sphere where 
the light cones point toward the singularity of the black 
hole, which automatically implies that material particles ap- 
proaching such boundary have no other fate than moving 
also toward the singularity and there is no need to impose 
a boundary condition there; this is a well known numeri- 



4 RESULTS 
Morphology 

We have found similar results as done in the past, par- 
ticularly, in the supersonic case, a shock cone character- 
ized by a zone of high density shows up at the rare part 
of the black hole. This behaviour appears to be a sort 
of late-time attractor behavior. In order to compare with 
previous research we show various examples. As illustra- 
tive cases we show in Figure [T] two examples of the shock 
cone morphology using KS coordinates, for two different 
values of the adiabatic index F = 4/3 and F = 5/3. 
Several other cases were studied by changing the veloc- 
ity of the wind and the rotation parameter of the black 
hole, see Table 1. The morphology of the rest mass den- 
sity is consistent with previous analyses (|Petrich et al. 19891 
IFont fc Ibafiez et al. 1998bl IFont et al. 1999|) . 

Another important check is the measure of the mass ac- 
cretion rate across a circle in the equatorial plane 9 = it/2, 
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Figure 1. Here we show the contour Unes of the rest mass density 
in order to compare the shock cone opening angle for models KS2 
(top) and KS8 (bottom) in Table. We then show that the angle of 
the cone is bigger for bigger values of the adiabatic index (Table) . 
The snapshot is shown for t = lOOOOM. The radial coordinate in 
units of the accretion radius race- 



M = — J^^ ay/jD{v'^ — P"^ /a)dcj>, as an indicator of the ac- 
cretion process during the evolution. We measure the mass 
accretion rate in terms of the proper time t at a given detec- 
tor. In Figure[2]we show the accretion rate for three different 
values of the initial velocity of the gas in order to illustrate 
the dependence of the accretion rate on the velocity of the 
wind. Our measurements are performed at a detector lo- 
cated at r = 2.1 for three different models. The figure shows 
that the higher the velocity of the wind the smaller the ac- 
cretion rate. Again these results are consistent with previ- 
ous studies (|Petrich et al. 19891 [Font fc Ibanez et al. 1998bl 
IFont et al.T999|l . 

Flip-flop behavior 

The flip-flop motion of the shock cone consists in the os- 
cillation of the cone itself along the angular direction. We 
illustrate in Fig. O snapshots of the shock cone morphology 
at various times for the model BL2. It can be clearly ob- 
served that the shock cone oscillates considerably in the an- 
gular direction. Similar results in the relativistic regime were 
reported recently for the flrst time in (|Donmez et al. 201ip . 

Something that called our attention is the possibility 
that the flip-flop behaviour might be due to the numeri- 
cal implementation of boundary conditions near the event 
horizon, because both in (|Donmez et al. 2011|l and in Fig.O 




Model KS1 
Model KS2 - 
Model KS3 - 



Figure 2. The mass accretion rates for models KS1,KS2,KS3. We 
show how the accretion process depends on the velocity of the gas, 
being the accretion rate higher for slower winds, as expected. 
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Figure 3. Snapshots of the rest mass density for the model BL2 
at various values of the coordinate time (top-left panel t = 7686, 
top-right panel t = 8062, bottom-left panel t = 8308 and bottom- 
right panel t = 8670) are shown. The motion of the shock cone 
along the angular coordinate is called the flip-flop behavior. We 
claim in this paper that the amplitude of such oscillations are 
highly dependent on the coordinates used to describe the black 
hole space-time, because they are singular at the event horizon. 
The units of the spatial coordinates are rescaled with race- 



BL singular coordinates at the event horizon are being used. 
Particularly important is that in order to simulate accretion 
processes it is necessary to use a radial domain such that 
G [rexc,rmax] with Tei^c > rEH, which implies that the 
inner boundary is time-like; moreover, if one wants to track 
the evolution of the gas up to regions near the black hole, 
the metric functions in non-penetrating coordinates tend to 
diverge and the gradients of the hydrodynamical variables 
are high. Then errors coming due to the implementation of 
a time-like artiflcial boundary near the event horizon are 
expected to contaminate the simulations in the numerical 
domain. 

In order to investigate this, we measure the oscillations 
along the angular direction of the shock cone when penetrat- 
ing KS coordinates are used. In Fig.|4]we show snapshots of 
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Figure 4. Snapshots of the rest mass density of the shock cone, 
using KS coordinates for the model KS2, are shown. The cor- 
responding times associated with each panel are: top-left panel 
t = 7686, top-right panel t = 8062, bottom-left panel t = 8308 
and bottom-right panel t = 8670. As we can see, the motion of 
the shock cone is very small compared with the case in which BL 
coordinates are used, see Fig. [S] By using the same coordinate 
time values as for the BL case to expose the snapshots, we do 
not intend to imply that the differences in the morphology are 
coordinate invariant. Instead, in order to compare the motion of 
the shock cones in the two coordinate systems we use a detector 
and the proper time as measured there as explained below. The 
units of the spatial coordinates arc rescalcd with race- 



what happens for the model KS2, physically equivalent to 
that in Fig. [3l The oscillations are rather small compared 
to those observed when BL coordinates are used. In Fig. [5] 
we show the amplitude of the oscillations when using the 
two coordinate systems for two representative cases in our 
Table. What we measure is the position of the maximum of 
the rest mass density along the angular coordinate in radi- 
ans on a circle of radius r = lOM; the maximum is located 
within the shock cone and is therefore a good quantity to be 
monitored; we plot it versus the proper time measured by a 
detector located at r = lOM, which appears in Figs. [3] and 
[4] as a circle. 

We identify flip-flop oscillations as those that clearly 
show the motion of the shock cone in the angular direc- 
tion. In Fig. [5] we distinguish at least two important modes 
when using BL coordinates, one with high frequency and 
small amplitude and the other with big amplitude and low 
frequency. The later is the one associated to the flip-flop mo- 
tion. The important observation here, is that when using KS 
coordinates, only the high frequency mode is observed, that 
is, the flip-flop mode does not happen. This clearly indicates 
that coordinates influence the motion of the shock cone, or 
what can be more interesting is that a true process of accre- 
tion, that is, allowing the gas to really enter the black hole, 
implies that the flip-flop effect does not take place. 

In order to check the accuracy of the numerical code, 
in Fig[6]we show the self-convergence of po using three reso- 
lutions for the case of KS coordinates. The self-convergence 
was calculated for a fixed angle within the shock cone. This 
is a very strong test because along such line the shock cone 
vibrates. We have a self-convergence better than first order 
when the shock cone has approached a nearly stationary 
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Figure 5. In these figures, we compare the position of the max- 
imum of the rest mass density along the angular coordinate <j) 
during the evolution in terms of the proper time. The models we 
are considering to carry out this comparison are: top panel: BL2 
vs KS2 and bottom panel BL5 vs KS5. In both cases the max- 
imum of Po is measured in a detector located at r = lOM. It 
is clear to see, from these figures, that when BL coordinates are 
considered, oscillations of high amplitude are presented, unlike in 
the case of KS coordinates. <j> coordinate is measured in radians. 



stage, while first order is acceptable for systems that include 
the evolution of shocks. 

In Fig. Owe show that using KS coordinates the oscilla- 
tions of the shock cone are not of the flip-flop type, instead 
they seem to be vibrations. We go further and look at the 
amplitude of such oscillations. In Fig.[7]we show the trend of 
the amplitude of the vibrations of the shock cone with reso- 
lution. In order to have an idea of what would happen in the 
continuum limit, we calculate the average of the amplitude 
of the vibrations for the two resolutions and measure the am- 
plitude of the vibrations with respect to such average. The 
result for the exercise in Fig [7] is that the amplitude around 
the average is 2.789e-2 for (Ar = 0.022, A0 = 0.0093) and 
1.08e-2 for resolution (Ar = 0.0146, A0 = 0.0062). 

In order to study the effects observed when using BL 
coordinates, we show in Fig.|8]the position of the maximum 
of the density within the shock cone in time along the an- 
gular coordinate at a fixed radius. We ask whether the am- 
plitude of the oscillations changes with resolution, and we 
find that the fiip-flop behaviour is independent of the reso- 
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Figure 6. Wc show the order of sclf-convergoncc Q of po for a 
fixed angle centered at the shock cone. Wc calculate the Li norm 
of the differences between the value of the density for the various 
resolutions. We show the self-convergence factor for resolutions 
(Ari,A0i) = (0.033,0.014), (Ar2,A</>2) = (An, A0i)/1.5, and 
(Ara, A03) = (An, A</,i)/(1.5)2. 
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Figure 8. We show the position of the maximum of the shock 
cone at r = lOM for various resolutions; we observe that the flip- 
flop effect persists, however the oscillations change the phase for 
the BL2 model. 
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Figure 7. The amplitude of the vibrations of the shock cone 
around the average and show how it decreases with resolution. 
We show only the oscillations of the maximum of the shock cone 
at a detector located at r = lOAf for the two highest resolu- 
tions we used. The oscillations start at around t ~ 500. Before 
that time we see a transient lapse that corresponds to the time it 
takes the shock cone to approach a nearly stationary regime. The 
top panel shows the result using medium resolution (Ar = 0.022, 
A(f) = 0.0093) and the bottom panel corresponds to high resolu- 
tion (Ar = 0.0146, A<f> = 0.0062). 

lution, as one may expect it could vanish in the continuum 
limit. The flip-flop behaviour is seen to be triggered with dif- 
ferent phases for different resolutions, which indicates that 
numerical discretization errors influence the appearance of 
the flip-flop which actually behaves as an instability. 



5 DISCUSSION AND CONCLUSIONS 

We study numerically the relativistic Bondi-Hoyle accretion 
on the equatorial plane on a space-time corresponding to a 
rotating black hole. Our analysis is focused on the shock cone 



motion of a supersonic gas that is being accreted onto the 
rotating black hole, and its dependence on the coordinates 
used to describe the space-time. We compare the cases of 
using BL non-penetrating horizon coordinates as done in 
previous research and KS penetrating coordinates. 

We flnd that when BL coordinates are used, the shock 
cone oscillates in a flip-flop fashion with amplitude of the or- 
der of a radian, whereas when using penetrating coordinates 
such behaviour does not take place. In fact, considering the 
frequency of the flip-flop oscillations when using BL coordi- 
nates, the oscillations observed when using KS coordinates 
do not seem to be of the same nature. Instead, in the case 
of KS coordinates we interpret the oscillations within the 
shock cone to be of vibration or resonant type. 

Furthermore, using various resolutions in our calcula- 
tions, we show that the amplitude of the shock cone vibra- 
tions when using KS coordinates decreases with resolution, 
suggesting thus that even the vibrations are not excited in 
the continuum limit. Our results are naturally restricted by 
the s-lab symmetry that we use in our calculations, that 
however has been also used to identify the shock cone oscil- 
lations with -for instance- QPOs (|Donmez et al. 201ip . We 
suggest in this paper that perhaps such oscillations are not 
physically sustainable, at least within s-lab symmetry and 
that 3D analyses are required. 

It is not possible to conclude about the condition that 
prevents the flip-flip effect from happening when using KS 
coordinates and excision inside the event horizon. We con- 
sider that allowing the matter to truly enter the black hole in 
the relativistic regime is the correct procedure. It may have 
important implications that have never been taken into ac- 
count before, specially in astrophysical scenarios, where usu- 
ally pseudopotentials of the Paczyhski-Wiita type are used 
to supplant a black hole, and usually it is assumed that mat- 
ter will be accreted when entering the accretion radius, and 
that furthermore require the implementation of an artificial 
time-like boundary and approximate boundary conditions. 
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